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Abstract 

The Homotopy Analysis Method (HAM) is employed in the analysis of non-linear one-dimensional heat 
conduction problem with temperature-dependent thermal conductivity. The temperature disuibution is obtained 
in terms of the scale space variable x and as a function of B which is a parameter in the heat conduction 
equation. It is observed that the parameter B has a strong influence over the rate of heat conduction in the bar. 
By choosing the convergence parameter, h, in a suitable way, we obtained solutions for the temperature 
distribution for various values of B. From this temperature distribution other heat transfer quantities can be 
obtained. It is observed that for some values of B the temperature distribution in the bar is an increasing function 
of B while for other values of B it a decreasing function. The results are displayed in figures. 
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Nomenclature 


A ; 


cross-sectional area of the material 


X 


dimensional space coordinate 


Dn, 


mth-order Homotopy-derivative 


K 


dimensional thermal conductivity 


k 


dimensionless thermal conductivity 


L 


length of the metal bar 


£ 


auxiliary linear operator 


T 


temperature 


A 


dimensionless space coordinate 




auxiliary parameter 


HC:J 


auxiliary function 


q 


embedding 


u 


dimensionless temperature 


B 


parameter in equation (3) 


Subscripts 


■3. 


ambient condition 


b 


end of the bar 


m 


mth order of approximation 


Greeks 




%m 


two- valued function 



q similarity variables 
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I. INTRODUCTION 

Nonlinear phenomena play a crucial role in applied mathematics, engineering and physics. The theory of 
nonlinear problems has recently undergone many studies. One of the nonlinear problems encountered in real life 
is the heat conduction in a bar with temperature dependent thermal conductivity. The thermal conductivity of a 
material generally varies with temperature. However, this variation is mild for many materials in the range of 
practical interest and can be disregarded. Even in such cases, an average value for the thermal conductivity is 
used and treated as constant. When the variation of thermal conductivity with temperature in a specified 
temperature interval is large; however, it becomes necessary to account for this variation to minimize error. 
Accounting for the variation of the thermal conductivity with temperature, in general complicates the analysis. 
Due to non-linearity of the problem, in most cases numerical approaches are common but analytical approaches 
seem to be intractable. The basic idea of Homotopy in topology provided an idea to propose a general analytic 
method for non-linear problems namely homotopy analysis method (HAM). The method was proposed by Liao 
SJ. [8] in 1992. This method is now widely used to solve many types of nonlinear problems. For example, 
Khani, F. et al, [7] used HAM for the solutions and efficiency of the non-linear fin problem in which he showed 
the efficiency of HAM. Abbasbandy, S. [3] applied HAM to solve a generalized Hirota-Satsuma coupled KdV 
equation. Liao et al [5] applied HAM for the solutions of the Blasius viscous flow problems. Other researchers 
have since employed HAM for solutions of seeming difficult scientific and engineering problems. In this paper, 
we introduce the basic idea of HAM and apply it to find an approximate analytical solution of the nonlinear 
one-dimensional heat conduction problem with temperature dependent thermal conductivity. 

II. MATHEMATICAL FORMULATION 

Here we consider a straight one-dimensional metal bar with a constant cross-sectional area. A,,. The bar 
with length, L, is exposed to a temperature source and extends into a fluid at ambient temperature T a . We 
assume that the conduction heat transfer is in a steady-state and no heat is generated within the bar. The one- 
dimensional steady state heat balance equation in dimensional form may be considered as 
d / dT, 

^feH' 0<X<L 

Where K is the thermal conductivity of the metal bar and in many engineering applications it may be assumed to 
be a linear function of temperature 

K = K(T] = K Jl + 1(1 - TJ] . f2) 
Here, K a is the thermal conductivity of the metal at the ambient temperature and X is the parameter describing 
the dependence on the temperature and is called the temperature coefficient of thermal conductivity. By 
introducing the following dimensionless variables 

u = i^t' x = i' k = Y 1 ' B=J ' (T t- T ^ & 

the general equation becomes 

-(k(U>-) = 0, 0<x<L, (4) 

Which is subject to 

U'(C0 =0 atx = 0 

U(l) =0 atx = L (5) 
The heat balance equation (4) can be rewritten as 

U"+BU"[J + BU 1 = 0 (6) 
In general, the analytic solution of Equation (6) is difficult. As a matter of fact, only if E = 0 will this equation 
posses an analytical solution. In engineering applications, an explicit approximation with an acceptable accuracy 
is preferred rather than an exact result in an implicit form with a rigorous proof. As a result, this paper intends to 
look for a solution algorithm which can generate an explicit approximate solution in terms of ordinary functions. 

III. A REVIEW OF HAM 

The basic idea of HAM which was introduced by Liao [8] is needed to deal with the goal of this work. The 
following definitions are needed for the start. 

Definition 1. Let <|) be a function of the Homotopy-parameter q. Then 
1 d m <k\ 

is called the mth-order Homotopy-derivative of <|), where m >0 is an integer. 
Let us consider the following nonlinear equation in a general form: 
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where 3f is a nonlinear operator, is an unknown function, and is temporal independent variables. For 
simplicity, we ignore all boundary or initial conditions, which can be treated in a similar way. 
This is a generalization of the traditional Homotopy method; Liao [8] in his paper constructs the so called zero- 
order deformation equation 

(1 - q)£[<Kg q) - u D (Tl ] = qhH <£) Mfyiq : q)I (8) 
where q E [0,1] is the embedding parameter, h is non-zero and is called the convergence parameter, H(£} is an 
auxiliary function, L is an auxiliary operator, Q[i(:J is an initial guess for and is an unknown function. A 
great freedom in choosing the auxiliary variables in Ham is of paramount importance. It is clear that when 
q = 0 and = 1 , we have 

0) = Uo® and+Gg; 1} = u(£>. 
Thus, as q increases from 0 to 1, the solution $(z: q] varies from the initial guess ^(z) to the solution utq}, 
Expanding $(z: q] in Taylor's series with respect to q, gives 

^:q)=u 0 © + ^u ln raq m (?) 
m=l 

where 

u m (fe) = D nl [^:q)]. (10) 
If the auxiliary linear operator, the initial guess, the auxiliary parameter, and the auxiliary function are properly 
chosen, the series (9) converges at q = 1 and we get, 

<l(z) = U D (y -f I^lUmCy. 
which must be one of the solutions of the original nonlinear equation, as proved by Liao[8]. For h = —land 
H(:J) = —1, equation (8) becomes 

(l - q M^ - u D <9] - q-W"l4><& q)l = 0- 

while the solution is obtained directly and without using Taylor's series. When l\(z) =1 Equation (8) changes 
to 

(1 - q)£l4(;; q] - Q D (y] = qbJf q)L (11) 
According to definition (10), the governing equation can be deduced from the zero-order deformation Equation 

(8), 

Define the vector 

5,1= {llfl (£), Ui(£) Un (4)}. 

Operating on both side of Equation (8) byD m , we have the so called mth-order deformation equation 

£tv ai 1&-Z ]B v a _ l QQ] = MK^R^^.Q (12) 

where 

Rju m _ v O = D m _ L (^(4; q)]), (13) 

and 

_ rO, m <1. 
m ll, m>2. 



Substituting (9) into (13), we have 

Rn-Gw-y = Dn,-i(MS; = E DUll ©q n ]l fl = D ) (14) 
It is of paramount importance to emphasize that u m for m > 1 is governed by the linear Equation (12) with 
the linear boundary conditions which comes from the original problem and can be easily solved by symbolic 
computation software, like Mathematical, Maple, and Matlab. 

IV. SOLUTION FOR THE TEMPERATURE DISTRIBUTION 

In this section, we obtain analytical approximate solutions of Equation (6), and to apply HAM to Equation 
(6), it is necessary to introduce the Molabahrami and Khani's Theorems. Molabahrami and Khani [6] proved the 
following theorems. 

Theorem 1. For Homotopy-series $ = H^=plJ n q a .. it holds that 

m 'i *i Pk_ 3 r^-j 

1-4=0 r 2 = D r 3 =D IV-2=0 lli-i=0 

wherem > 0 and k = 0 are positive integers. 
Proof: See [6] 



www.theijes.com 



The IJES 



Page 39 



Homotopy Analysis Method for One. 



Corollary l.from theorem 1 , we have 

m r i r 2 rit-3 rk-i 

Dm^ -1 *') = ^ w m _ ri ^ w ri _ r: ^ w rj _ l3 ... ^ w^^^ "Wn^w 
r u =0 r : = D r 3 = P Tk-j=B rjj-^O 

and 

m r i r 2 Ht-s rk-! 

r t = D r2=B r 3 = B rk_2 = D rt_ 1 = P 

from initial conditions (5), it is reasonable to express the solution by the set of base functions 

{enCxb > 0)1 (15) 
in the form 

U(x) = ^ a^Cx), (16) 

□ = D 

wherea D is a coefficient and e n (x) is determined to provide the so-called rule of solution of expression u(x) . 
Under the rule of solution denoted by (16), and from Equation (6), it is straightforward to choose 

Ll$(x:q)]= (17) 
ox- 1 

as an auxiliary linear operator which has the property 
£[c 0 + c L x] = 0 

where c p and Cj_ are coefficients. The solution given by HAM denoted by (16) can be represented by many 
different base functions. In this work, we set e n (x) = x :n . According to Eq. (6), we define the non-linear 
operator 

M*(x; q]] = B^Cx; q) +■ D ( ^f-) (13) 

From Equations (13) and (18) and Theorem 1 and Corollary 1, we have 

= Dm-iC^BD^M") + B D m _ ]_((<(>} 
From the above equation, the following results are achieved: 
KiFM = U B +BW + B(U D y, 
R s [l? L (x)] = U L " + 5(0,0," + U L U D ") + B(2U 0 'C l l ), 
RjPzCjO] = U 2 " +B(U D U 2 " + U 1 U 1 " + U 2 U D ")-KB(2U D r U/+(U 1 ^ 2 ). 
etc., where the primes denote differentiation with respect to x. 

Now, the solution of the mth-order deformation Equation (11) under Equation (18) with initial conditions 
U' m (0) = 0 and U m (l) = 0, for m > 1, is 

U m W = Z m U m _ L (£ 4- h / D X dw / D W H(v)R m (^(v)) dv + c lX + c D . (19) 

According to Equations. (6) and (5) and the rule of solution expression (16), it is straightforward to show that 
the initial approximation can be written in the form 

U fl (x)=x. (20) 
According to the rule of expression denoted by (16) and from Equation (19) the auxiliary function H (x) , should 
be as follows; 

HOO = l 

Now from Equation (19), we can successively obtain 
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3 

h Bh 1 2 hx ' 

U , (x ) = - — - + — Bh x + 

6 2 2 6 

h Bh h 2 2 5/i 2 5 2 A 2 1 , , hx 3 1 , , , 1 , , 

[/ (*) = - — - - - - + —Bh'x' + + —B h x + — Bh x 

6263 22 62 6 



U 3 (x)- 



h Bh h 2 A Bh 2 , , h 3 53 S/i 3 41 B "V B 3 /z 3 
— - B h - — 



+ — Bhx 2 + Bh 2 x 2 + — 



3 7 

Bh x 



72 



60 



1 

12 



2 3 2 1 3 3 2 h X 2 2 3 

5 jc - — B h x + + + B h x ' + 



k 3 X 3 1 33 7 2 ,3 1 

Bh x + — B h x + — 



Bh 3 x 4 



6 36 



12 



11 



60 



d 2 ; 3 5 
Bhx 



nr 3 6 

+ Bh x 

72 



In this way, we derive U m (x] for m =1,2,3 successively .At the k th -order approximation, we have the 

analytic solution of 
k 

V(x) * y LT m U). 



analytic solution of Equation (6), namely 



(21] 



Liao [8] presented the auxiliary parameter in order to control the rate of convergence of the approximate 
solutions obtained by HAM. The parameter h. determines the convergence region and the rate of approximation 
for HAM. For this purpose, h. — curves are displayed in Figures 1-3. It is clear from these figures that for values 
of B < 0 and B >l,h= —1 fall within the convergence region, i.e., that HAM is convergent for these values 
of B. Considering figures 1-3, it is further observed that h converges at (-1.8 < h < -0.4), therefore taking h= - 
0.5 enables the plotting of U(x) against x for various values of B in figures 4-12. 

V. RESULTS AND DISCUSSION 

From figure 4, it seen that the temperature distribution when (B=0) is a linear function of x. This is the 
simplest expression of the temperature distribution in the bar. 

Figures 5, 6, 8 and 12 illustrate the temperature distribution for various values of B. It is clear from the 
figures that the temperature distribution in the bar is an increasing function of B, while in figures 7, 9, 
10 and 1 1 the temperature distribution in the bar is a decreasing function of B. 
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Figure 1. The h-curves for B=0.4, for 4 th -order approximation of U" [0] and U 4 [0]. 
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Figure 5. Temperature distribution for B= 0.1, 0.2, 0.3. 
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Figure 6. Temperature distribution for B= -0.1, -0.2, - 0.3 
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Figure 7. Temperature distribution for B= 0.4, 0.5, 1. 
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Figure 8. Temperature distribution for B= -0.4, -0.5, -1. 
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Figure 9. Temperature distribution for B= 0.7, 0.8, 0.9. 
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Figure 10. Temperature distribution for B= -0.7, -0.8, -0.9. 
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Figure 11. Temperature distribution when B= 1.2, 1.5, 2. 
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Figure 12. Temperature distribution for B= -1.2, -1.5, -2. 
VI. CONCLUSION 

This work employs HAM to solve the nonlinear steady state heat conduction problem with temperature 
dependent thermal conductivity. The results show that B has a very strong influence on the temperature 
distribution in the bar and cannot always be ignored without error. For some values of B, the temperature 
distribution in the bar is an increasing function of B while for other values of B it a decreasing function. The 
results are display in figures. From the analytic expression in equation (21), other heat transfer quantities can be 
obtained easily. 
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